(N 



> 



X 



Precise Neutron Magnetic Form Factors 



-k 



G. Kubon^ H. Anklin^ P. Bartsch*^ D. Baumann*^ 
W.U. Boeglin^'i K. Bohinc^ R. Bohm^ M.O. Distler^ 
I. Ewald'^ J. Friedrich*^ J.M Friedrich*^ M. Hauger^ 
^ ■ A. Honegger^ P. Jennewein"^ J. Jourdan*^'^ M. Kahrau*^ 

(N '• K.W. Krygier^ A. Liesenfeld^ H. MerkeP U. Miiller^ 

^ ; R. Neuhausen*^ Ch. Normand^ Th. Petitjean^ Th. Pospischil^ 

;z; : M. Potokar^ D. Rohe^ G. Rosner"'^ H. Schmieden^ I. Sick^ 

oo ■ S. Sirca^ Ph. Trueb^ A. Wagner*^ Th. Walcher*^ G. Warren '^ 

M. Weis^ H. Wohrle^ M. Zeier^'^ J. Zhao^ B. Zihlmann^'^ 



^Dept. fiir Physik und Astronomie, Universitdt Basel, Basel, Switzerland 
^Institute Jozef Stefan, University of Ljubljana, Ljubljana, Slovenia 



o 

C^^ ■ ^Insitut fiir Kernphysik, Johannes Gutenberg-Universitdt, Mainz, Germany 

O 

o 

^- 

QJ , Abstract 

^ , Precise data on the neutron magnetic form factor Gmn have been obtained with 

rH I measurements of the ratio of cross sections of D(e,e'n) and D(e,e'p) up to mo- 

mentum transfers of Q^= 0.9 (GeV/c) . Data with typical uncertainties of 1.5% 
are presented. These data allow for the first time to extract a precise value of the 
magnetic radius of the neutron. 
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Introduction: Detailed information on the inner structure of the nucleon is 
provided by accurate data on the dependence of the nucleon form factors on 
momentum transfer Q^. Such data serve as a sensitive test for models of the 
nucleon. Particularly at low Q^ accurate data on the form factors allow for 
both a determination of the electromagnetic radii and accurate calculations 
of nuclear form factors. 

While the proton form factors are known with excellent precision over a large 
range of Q^, data for the neutron are of much poorer quality due to the lack of 
a free neutron target. This is true for both the electric form factor, Gen, and 
to a somewhat lesser extent for the magnetic one, Gmn- However, with today's 
high-duty factor, high-current electron beam facilities and the advances made 
in polarized beam and target technology, significant progress in this area is 
being made. 

In the past, Gmn has been determined mostly from quasi-elastic D{e, e') cross 
sections (see references in [1]). The extraction of Gmn requires a longitudi- 
nal/transverse separation and a subtraction of the (dominant) proton mag- 
netic contribution. The uncertainties resulting from the deuteron wave func- 
tion, meson exchange currents (MEC), and final state interactions (FSI) are 
greatly amplified by the two subsequent subtractions and limit the accuracy 
of Gmn from these experiments to ~20%. Due to these limitations alternative 
techniques have been used in recent experiments. 

One of the techniques determines Gmn from the asymmetry measured in 
^if e(e, e')-scattering [2,3]. This challenging technique is presently limited to 
low Q^ where today's rigorous non-relativistic 3-body calculations can be 
applied to remove the dependence on the nuclear structure, FSI, and MEC 

[4]. 

The neutron magnetic form factor can also be obtained from an exclusive 
cross section measurement of D{e,e'n) [5]. This technique avoids essentially 
the subtraction of the proton contribution which was responsible for a part 
of the large sensitivity to systematic errors in the past. However, the method 
still depends on a deuteron model for the extraction of Gmn- 

The best method to minimize the sensitivity to the nuclear structure is a 
determination from the ratio R = (7{e, e'n)/a{e, e'p) on the deuteron in quasi- 
free kinematics [1,6-8]. The ratio is insensitive to the deuteron wave function 
and corrections due to FSI and MEC are calculable and small. The price 
to pay is the need for a precise measurement of the absolute efficiency rj of 
the neutron detector employed. A measurement of rj and a detailed study of 
the detector response, however, is possible when using high-intensity neutron 
beams available at the proton-beam facilities [1,7]. 

In a pilot experiment it was demonstrated that this method leads to a deter- 



mination of Gmn with an accuracy of 1.7% [7]. Similar measurements over an 
extended Q^-region are possible with a high-duty factor electron accelerator 
like the Mainz Microtron (MAMI) [9]. Precise measurements of the ratio R 
were performed in the Q^-range from 0.2 to 0.6 (GeV/c) and values of Gmn 
were extracted with error bars as low as 1% [1]. In the present work we extend 
the Q2-range and present data of Gmn for Q^ of 0.071, 0.125, 0.359, 0.894 
(GeV/c) in the following labeled as kinematics 1 to 4. To check the consis- 
tency of such measurements the point of the pilot experiment [7] has been 
re-measured (label 2) using a different nucleon detector, electron and neutron 
beam facility. 

Measurements at MAMI: At MAMI the yield of D{e, e'n) and D{e, e'p) in 
quasi-elastic kinematics was measured with electron beam energies of 600 MeV 
and 555 MeV for kinematics 1 and 2 and of 855 MeV for kinematics 3 and 4. A 
beam current of 0.5 fiA incident on a cylindrical 2 cm thick liquid deuterium 
target cell with 7 ^m HAVAR windows was employed for kinematics 2 to 4. 
Due to the low proton energy (36 MeV) a target cell with a lateral width of only 
1 cm liquid deuterium was used for kinematics 1 in order to minimize multiple 
scattering and energy loss effects of the knocked-out protons. Spectrometer A 
[10] with a solid angle of 28 msr detected the scattered electron in coincidence 
with the recoiling nucleon. 

The nucleon detector consisted of two 10 cm thick plastic converters, Ef and 
Efy, preceeded by 3 thin AE counters used to identify the incident nucleon. 
The thickness of the AE counters was 1.5 mm (5 mm) for kinematics 1 to 2 
(3 to 4). Except for the opening towards the target, the detector was shielded 
with 10 cm thick lead walls. Lead absorbers of 0,1,3, and 20 mm thickness 
were placed at the entrance window of the detector for kinematics 1 to 4. The 
absorbers are needed in the D{e,e'n) measurements in order to absorb low 
energy photons. The nucleon detector covered an angular range of ±78 mr in 
horizontal as well as vertical direction resulting in a solid angle of 24.3 msr. 

The yield ratio has been determined via the simultaneous measurement of the 
D{e, e'n) and the D{e, e'p) reactions which makes the ratio independent of the 
luminosity, dead time effects, and the efficiency of the electron arm. Even at 
the highest Q^ a signal-to-noise ratio of > 50 for the (e, e'n) measurement 
has been achieved. 

In the analysis neutrons were defined requiring no hit in at least 2 of the 3 
veto counters {AEi ■ AEj with i,j = 1, 2, 3) and a hit in Ef in the coincident 
time window. The threshold used in the definition of neutrons in Ef was set at 
25 MeV (10 MeV) for kinematics 2 to 4 (1). The number of counted neutrons 
was corrected for the efficiency of the veto condition and for misidentified 
protons due to inefficiencies of the veto counters. The veto efficiency correction 
was determined to <1% for kinematics 1 and 2, 11.1% (30.4%) for kinematics 



3 (4). The correction for misidentified protons ranged from 0.8% to 5.2%. The 
three independent measurements for the number of neutrons, determined with 
the three different veto conditions, agreed within 0.16%(1.5%) for kinematics 
1 to 3 (4). The agreement of the neutron counts within 1.5% for kinematics 4 
provides confidence on the vahdity of the rather large veto correction. 

Protons were counted in the TOF spectrum of Ef with a coincident time 
signal in at least one of the three AE counters. The number of protons was 
corrected for the inefficiency of the AE counter used. The final number of 
protons, obtained using different AE,Ef combinations, agreed within 0.11% 
for kinematics 1 to 3 and to 0.2% for kinematics 4. The uncertainties in the 
corrections to the number of protons and neutrons are small relative to the 
statistical error and are included in the error of the yield ratios which are 
summarized in table 1. 

Measurements at PSI: The determination of the neutron detection effi- 
ciency was performed using the high energy neutron beam (100-590 MeV) at 
the Paul Scherrer Institut (PSI) [11] for kinematics 2 to 4 and the monoener- 
getic neutron beam of 68 MeV for kinematics 1 [12]. 

The high (low) energy neutron beam is produced via the C(p,n) (D(p,n)) - 
reaction with a neutron flux of 10^/s (10^/s) for a 5 ^A proton beam. A tagged 
high intensity neutron beam was produced via the H{n,p)n reaction, scatter- 
ing the neutrons from a 1 cm thick liquid H2-taiget. The recoil protons were 
detected with AEp and Ep plastic scintillators recording the amplitude and 
the time-of-flight (TOF). Four multi-wire proportional chambers (MWPC) 
determined the proton trajectory, thereby flxing the target coordinates with 
an accuracy of ±1.5 mm and the recoil angles with an accuracy of ±0.2°. 

Together with the measured time reference of the radio frequency of the cy- 
clotron, the recorded information allowed for a determination of the incident 
neutron energy, the energy of the recoiling proton and its recoil angle. This 
provided a redundant determination of the energy and position of the tagged 
neutron beam which was free from contributions of background reactions. 

The accuracy of the tagged neutron energy from 100 to 490 MeV ranged from 
1.9 to 16.6 MeV. For 36 MeV tagged neutrons the energy was determined to 
an accuracy of 0.2 MeV. The position on the detector surface was determined 
with an accuracy of ±3.5 mm. 

The nucleon detector was placed in the tagged neutron beam, and its absolute 
efficiency distribution ri{x,y,Tn) was measured as a function of the point of 
impact {x,y) and the neutron energy (T„). The knowledge of the distribu- 
tion ri{x,y,Tn) is necessary because the detector illumination is unavoidably 
different for the H{n,p)n and the D{e,e'n) data taking. Use of the identical 
configuration during the ratio measurements and the efficiency measurements 



ensured that absorption effects of neutrons were automatically included in the 
efficiency determinations. A more detailed account on the neutron efficiency 
determinations is given in [13,14,6]. 

Considerable care was taken to ensure identical cuts on the energy deposited 
in the detectors. H{e,e'p) at MAMI and H{n,p) at PSI were used to obtain 
absolute calibrations of the amplitude spectrum. In addition, each converter 
was monitored for both gain variations and baseline shifts using two temper- 
ature compensated LED's whose light output was in turn monitored by very 
stable PIN-diodes [7]. This information allowed for a determination of the en- 
ergy scale of the converters to an accuracy of 0.8% during the PSI and the 
MAMI runs. 

Further corrections, relevant in the determination of R were measured with 
tagged protons by placing the detector on the recoil arm. The dominant cor- 
rection of up to (12.9±0.7)% was due to multiple scattering and energy loss 
effects leading to proton losses in the lead absorber. In the same arrangement 
the variation of the light collection efficiency, required to match the measured 
ri{x,y,Tn) to the neutron D(e, e'n)-distribution, was measured. 

Two efficiency measurements bracketed the measurement of the yield ratio 
at MAMI in order to check the reliability of 77. Consistent results for 77 were 
found for all kinematic and absorber conditions. 

Evaluation of R: From the measured yield ratios and efficiencies the D{e, e'n) 
to -D(e, e'p) cross section ratios R^^^ were determined with accuracies of 1.7% 
to 2.7%. The resulting values of R^^^ are independent of the applied Ef thresh- 
old. The final results for R'^^ and the corrections applied in its determination, 
all of which are based on measurements, are summarized in table 1. 

To extract Gmn from the measured _R^^p values we have to take into account 
the effects of FSI, MEC, and isobar currents (IC) beyond the plane wave im- 
pulse approximation (PWIA). To this end, we used theoretical results oi R^^'^° 
calculated with the Paris potential by Arenhovel [15]. The total corrections are 
listed in table 2 as a deviation D = {Rp\yia/ R^^'^°) — 1 from the PWIA-value. 
The dominant contribution to D (up to 99%) is due to FSI, mostly charge 
exchange scattering. The contributions from MEC and IC are of order 0.5%, 
and relativistic effects are negligible [16]. 

At Q^ = 0.071 (GeV/c) where the correction is largest the dependence of D 
on the nucleon-nucleon potential was studied. A variation of typically 0.9% in 
D was found when using the Bonn-R-space, the Argonne V14, and the Paris 
potential. This is well within the relative systematic uncertainties (FSI: ±8%, 
MEC: ±40%, IC: ±60%) used to calculate the errors of table 2. 

The resulting value R^pwia = {o'e-n/o'e-p)pwiA = R^^^ ' (1 ± -D) is the exper- 



Q2 (GeV/c)2 


0.071 


0.125 


0.359 


0.894 


Yield ratio 


±0.5 


±1.0 


±0.9 


±2.1 


Illumin. matching 


±1.0 


±0.5 


±0.7 


±0.6 


Thresh, cahbration 


±0.8 


±0.8 


±0.8 


±0.8 


Proton losses 


-7.6±0.9 


-11.7±1.0 


-8.7±0.7 


-12.1±0.7 


(p,n)-correction 


-0.00±0.02 


-0.07±0.02 


-0.3±0.1 


-1.9±0.6 


Error of r] 


±0.6 


±1.0 


±0.5 


±0.8 


Tn and x, y 


±0.4 


±1.1 


±0.6 


±0.8 


R 


0.0973 


0.139 


0.279 


0.311 


Rel. error of i? in % 


±1.7 


±2.0 


±1.7 


±2.7 



Table 1 

Results for R and the experimentally determined errors and correc- 
tions in % of R. Except for the yield ratio and rj the errors are mainly 
systematic in nature. 

iinental ratio of the e — n cross section (which is essentially given by Gmn) to 
the e — p cross section corrected for non-PWIA contributions. The contribu- 
tion of the neutron electric form factor to R'pwia is small; it introduces only a 
small additional uncertainty despite the poor knowledge of Gen (see table 2). 



Q2 (GeV/c)2 


0.071 


0.125 


0.359 


0.894 


D{%) 


-24.1 ±2.0 


-9.9 ±0.9 


-3.7 ±0.4 


-1.2 ±0.2 


Contr. of Gen (%) 


-1.4 ±0.3 


-1.9 ±0.4 


-2.8 ±0.6 


-0.6 ±0.1 


ae-p/(JD 


0.948 ± 0.009 


0.931 ±0.014 


0.933 ±0.019 


1.089 ±0.023 


Gmn/ifJ-nGD) 


0.990 ±0.013 


0.967 ±0.013 


0.989 ± 0.014 


1.062 ±0.017 



Table 2 

Results for Gmn relative to the dipole form factor Go = (1 ±(5^/0.710)~^. 
The error on Gmn includes both the experimental contribution (table 1, 
(Te_p, Gen) and the ones due to theory. 

To evaluate the e — p cross section, o"e_p, we used the world's supply of ae-p 
data in a range of 0.5 fm~^ around the desired Q. In this range of Q, we used 
a parameterization with a relative Q-dependence as given by the Mergell et 
al. fit [17] to Gep and Gmp, with the overall normalizations fitted to the world 



data. The statistical errors of the data have been treated in the standard way, 
while the systematic errors have been accounted for by changing the data of 
each set by its error, refitting, and adding the changes due to systematic errors 
in quadrature. The resulting e —p cross sections relative to the ones obtained 
with dipole form factors are listed in table 2 together with the final results for 
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Fig. 1. The figure shows the present results (t), together with results of 
the previous MAMI/PSI [1] (a), the NIKHEF/PSI [7] (♦), the Bates ^He [2] 
(x), and the JLab ^He experiments [3] (•), in comparison to various model 
calculations. Solid: Mergell tt al. [17], dot: Kubis [18], dash-dot: Eich [19], 
dash: Schlumpf [20], dash-dot-dot: Lu et al. [21]. 

The present data on Gmn are shown in figure 1 together with the recently 
determined data sets (> 1990) of refs [1-3,7] and some recent calculations. The 
present data extend the previously investigated Q^-region in both directions 
and allow for a direct comparison of the data measured at NIKHEF/PSI [7] 
and JLab [3] . The agreement of these data measured with different techniques 
at different facilities at an unprecedented level of precision is very satisfactory. 

On the other hand, these combined data do not agree with the measurements 
by [5,8] (not shown in figure 1). This is due to the fact that in these ex- 
periments the three-body reactions D{e,p)ne' and H{e,TT)ne' were used to 
tag the recoiling neutron in the r^-determination. This, however, would have 
required significant corrections for neutrons that miss the nucleon detector; 
these corrections were not applied [22]. 

The data shown in figure 1 clearly differ from the crude empirical expression 
Gd = (1 + Q"^ /m^)~'^ with irP' = 0.71 (GeV/c) used to remove the dominant 
Q^-dependence in table 2 and figure 1. In addition, the data show significant 
differences to both the non-relativistic constituent quark model calculation by 
Eich [19], and the relativistic version by Schlumpf [20]. Similar differences are 



observed when comparing the data to a recent cloudy bag model calculation 
by Lu et al. [21]. The data are also compared to results of the relativistic 
chiral perturbation theory by Kubis and Meissner [18] and the calculation by 
Mergell et al. [17] based on a fit of the proton data using dispersion theoretical 
arguments. While none of these calculations describe the data satisfactorily 
the tendency is given by the calculation of Ref. [17]. 

Parameterization of Gmn- The present data allow for the first time a purely 
experimental extraction of the root-mean-square (rms) magnetic radius of the 
neutron defined as: 

< r^ >= 

^ mn ^ 




In the past, the experimental information on r^n was based on a dispersion- 
theoretical analysis of the combined set of electromagnetic form factors of 
both neutron and proton. In this framework Tmn is determined mostly from 
constraints other than the experimental Gmn data which up to now were very 
limited in accuracy. 

The present determination of Tmn uses only the experimental data on G„m- 
The data of figure 1 together with the higher Q^ data from [23] are taken into 
account. A continued fraction parameterization given by: 

1 + ^l^ 



1 + -- 



is fitted to the data. The magnetic radius of the neutron is related to the 
parameter bi via < r^^ >= 6 ■ bi. 

The continued fraction representation [24] has the advantage to converge in a 
wider domain than the usual power series expansion. One finds, for example, 
that five terms are sufficient to reproduce Gmn of the dispersion analysis by 
Hoehler [25] or Mergell [17] for Q^ up to 4 (GeV/c)^ 

Figure 2 shows the result of the fit using the five terms bi, ..65 = 3.26, —0.272, 
0.0123, —2.52, 2.55 (GeV/c)~ with a Xred = 0.91. For the magnetic radius 
we find Tmn = 0.873 ± 0.011 fm. The uncertainty of 1.1%, an improvement of 
a factor of 10 compared to previous determinations, covers the error of the fit 
(0.6%) and a systematic uncertainty due to the choice of the fit function. The 
latter was determined by fitting with different parameterizations pseudo data 
from a dispersion analysis placed at the Q^-values and with the errors of the 
real data. 
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Fig. 2. The figure shows the continued fraction fit to the data. Symbols for 
the data as in figure 1) plus the data by Lung et al. (+) [23]. 

The present result provides an additional observable to test theoretical cal- 
culations of the nucleon. The calculation by Buchmann et al. [26] has been 
performed in the framework of a constituent quark model including gluon, 
meson, and confinement exchange currents. The calculation predicts a radius 
of 0.891 fm in marginal agreement with the present experimental result. The 
agreement is worse when comparing to the very recent result of 0.84 fm by 
Kubis and Meissner [18] obtained in fourth order relativistic baryon chiral 
perturbation theory. 

Summary: The present data combined with the other data of figure 1 have 
improved our knowledge of the neutron magnetic form factor in the Q^-range 
from 0.07 to 0.89 (GeV/c) by a factor of 10 compared to determinations based 
on quasi-elastic (e, e')-data performed in the past [1]. The improvement is 
mostly due to the fact that the exploitation of the ratio R = a{e, e'n)/a{e, e'p) 
depends least on the input of theory; the measurement of R becomes possible 
by performing the needed calibrations using a high-intensity tagged neutron 
beam, and by using an electron beam with 100% duty factor. The increase of 
accuracy allows for a detailed comparison to theoretical model calculations. 
The agreement of three high-precision experiments gives us confidence in the 
reliability of the gained knowledge on the magnetic properties of the neutron. 
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